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Abstract. The cloud of cold atoms obtained from a magneto-optical trap is known to exhibit two types of 
instabilities in the regime of high atomic densities: stochastic instabilities and deterministic instabilities. 
In the present paper, the experimentally observed deterministic dynamics is described extensively. Three 
different behaviors are distinguished. All are cyclic, but not necessarily periodic. Indeed, some instabilities 
exhibit a cyclic behavior with an erratic return time. A one-dimensional stochastic model taking into ac- 
count the shadow effect is shown to be able to reproduce the experimental behavior, linking the instabilities 
to a several bifurcations. Erraticity of some of the regimes is shown to be induced by noise. 



PACS. 32.80.Pj Optical cooling of atoms; trapping 
ical systems - 05.40.Ca Noise 
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in two groups: stochastic and deterministic behaviors. For 
the former, the dynamics originate in noise, i.e. in dynam- 
ical components, usually with a large number of degrees 
of freedom, considered as external to the system. This is 
usually experimental technical noise, and requires to add 



1 Introduction 

Experimental quantum physics knows since several years 
spectacular results, thanks to a simplification to produce 
quantum objects with long coherence times or macroscopic 
dimensions. Let us cite the achievement of the Bose-Einstein in the model a random component. Such a complex dy- 
condensatespQ, the characterization of quantum chaos 0, namics is meaningless from the physical point of view, 
the improvement of atomic clocks the designs of quan- because it cannot give any new informations about the 
turn computers and quantum communication systems 0], MOT mechanisms. On the contrary, deterministic dynam- 
or also the accurate understanding of quantum decoherencej^cs are intrinsic to the system, and do not require to add 
One of the basic tools used to obtain most of these re- anything to the model: periodic instabilities can appear 



suits is the Magneto-Optical Trap (MOT), which performs 
the cooling of atoms at temperatures of the order of the 
/zK: this is the first step before reaching lower tempera- 
tures where the quantum properties of atoms dominate. 
Although it is a key device in the new atomic physics, 
the basic mechanisms determining the properties of the 
cloud of cold atoms in a MOT have been poorly studied, 
and the collective dynamics of these still "classical" atoms 
have been almost ignored, even though the existence of 
instabilities in the MOT is known since the first realiza- 
tions. On the contrary, some simple empirical rules are 
used to avoid these inconveniences. Nevertheless, an accu- 
rate knowledge of the individual and collective behaviors 
of the cold atoms in the cloud could help in understanding 
the limitations of the process, and above all, to enhance 
it through the control of the dynamics, as it was done in 
many other systems, in physics and other fields of science. 

However, a necessary preamble to such applications is 
the identification of the nature of the dynamics observed 
in MOTs. Indeed, complex behaviors may be subdivided 



with two degrees of freedom, while chaos needs at least 
three degrees of freedom. This last case opens many per- 
spectives: for example, it is possible to reach new working 
points by the methods of control of chaos, or to measure 
parameters which are usually inaccessible 



Recent studies have shown that the collective behav- 
ior of the atomic cloud produced by a MOT exhibit both 
stochastic instabilities^!!] and deterministic instabilities 0. 
The former has been extensively described in [5]. A model 
demonstrates that the different stochastic behaviors ob- 
served in the experiments are well explained if the absorp- 
tion of light by the atoms is taken into account, through 
the so-called shadow effect [2| ■ It is also shown that these 
stochastic instabilities are not "instabilities" in the usual 
meaning, as they result from an amplification of noise, 
due, from a dynamical point of view, to the folded struc- 
ture of the stationary solutions. The same model was also 
predicting, for slighly different values of the parameters, 
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deterministic instabilities which, in turn, has been ob- 
served experimentally 

In the present paper, we report an extensive study of 
these deterministic instabilities. We detail and complete 
the experimental results given in 2], and analyze accu- 
rately the mechanisms leading to the different determinis- 
tic regimes through the model introduced in [jy . We show 
in particular that the model is able to reproduce each 
type of dynamics, and predicts deterministic chaos. We 
show also the main role that noise is still playing in the 
dynamics. 

The paper is organized as follows. After this introduc- 
tion, section describes briefly the experimental setup, 
and section |3 gives a detailed analysis of the experimental 
observations. Section 01 is devoted to a short description 
of the model, already detailed in 9 . In section the sta- 
tionary solutions of the model are discussed, and in section 
El the deterministic dynamical behavior predicted by the 
model is described, and compared with the experiment 
results. Finally, in section [T] the effect of noise on the dy- 
namics is studied. 



2 Experimental set-up 

The experimental set-up has already been described in 
detail elsewhere [H], and thus the description here is sim- 
plified. The Cesium-atom MOT is in the usual cr+ — cr_ 
configuration, with three arms of two counter-propagating 
beams obtained from the same laser diode. The waist 
wt of the trap beams may be varied from typically 3 to 
10 mm. We use a configuration where counter-propagating 
beams result from the reflection of the three forward beams. 
This simplifies the detection of the dynamics, as compared 
to a six independent beams configuration. Indeed, because 
of the shadow effect, a center-of-mass motion is generated. 
However, as the nonlinearities involved in both cases are 
the same, we expect that the dynamics will be fundamen- 
tally of the same nature in the two configurations. 

The dynamics of the atomic cloud consists in a defor- 
mation of the spatial atomic distribution, leading to fluc- 
tuations of the shape of the cloud, as illustrated in Fig. 
^ Therefore, the relevant dynamical variables allowing us 
to describe instabilities, could be the shape of the cloud 
(i.e. for example the local velocities and atomic densities 
in the cloud). This type of description corresponds to a 
high dimensional model, associated with partial differen- 
tial equations. Here, for the sake of simplicity, we choose 
to limit our description to the center of mass (CM) lo- 
cation r, and the total number of atoms n in the atomic 
cloud. This allows us to model the system with ordinary 
differential equations, and reduces the dimension to seven, 
and even three in a ID model. As it is shown in the fol- 
lowing, the use of this description appears to be sufficient 
to understand the main mechanisms of the instabilities. 

A 4-quadrant photodiode (4QP) is used to detect the 
fluorescence of the cloud. The differential signal of the 4QP 
allows us to measure the motion of the CM through one of 
its component r, while the total signal gives us the number 
n of atoms inside the cloud. A second 4QP, perpendicular 



Fig. 1. Sequence of snapshots showing the time evolution 
of the unstable atomic cloud. Snapshots are presented in the 
chronological order, each one being separated by 40 ms. This 
sequence corresponds to the fast stage of aCp cycle (see sec- 
tion ej. 

to the first one, prevents the measure from line-of-sight ef- 
fects due to the optical thickness of the cloud. We checked 
that whatever the type of dynamical behavior, the signal 
recorded by both 4QP have the same properties and are 
qualitatively identical. 

Parameters acting on the dynamics have been exten- 
sively discussed in [Hj. The detuning A of the MOT, the 
magnetic field gradient G, the MOT beam intensity 1\ 
and the repumper laser intensity I rep may be considered 
as control parameters, because they can be easily changed 
in the experiment. On the contrary, the alignment of the 
MOT beams, the vapor pressure in the cell and the MOT 
beam waist, which also play a crucial role in the dynamics, 
cannot be considered as control parameters, either because 
they cannot be changed easily, or because they cannot 
be measured with accuracy. Therefore, these parameters 
have not been varied in the experiments. The parameter 
ranges explored in the present experiment are summarized 
in Tab. □ 

Table 1. Range of the parameters used in the present exper- 
iment. G is the magnetic field gradient, 7+ is the intensity of 
the forward beam and 6 is the detuning. I s is the saturation 
intensity ( I 3 — 1.1 mW) and r is the natural width of the tran- 
sition. The last column indicates the default parameter values 
used to obtain the results reported in the present paper. 

range default set 

G G < 14 Gem" 1 14 Gcnr 1 
7+ = I/I s 4 < 7+ < 20 10 
A = 5/r A<-0.5 



3 Experimental results 

In [J]], it has been shown that the atomic cloud exhibits two 
types of instabilities, depending on the parameters of the 
MOT, in particular the trap laser beam intensity I\. When 
Ji is small, typically less than IO/5 (I3 = 1.1 mW/cm 2 is 
the saturation intensity), instabilities are essentially stochas- 
tic. Depending on the other parameters, as the detuning, 
several types of stochastic instabilities occur. In the Sl 
behavior, instabilities are characterized by a unique slow 
time scale, and no component larger than 2 Hz appears 
neither in the motion of the trap, nor in its population. 
On the contrary, in the Sh behavior, a second time scale, 
at higher frequency (typically from 20 to 100 Hz) appear 
in the trap motion, but not in the population^ . 



Andrea di Stefano, Philippe Verkerk and Daniel Hennequin: Deterministic instabilities in the magneto-optical trap 



3 



When I\ is increased, S instabilities are progressively 
replaced by C instabilities ( C stands for cyclic) . These in- 
stabilities have already been described in |S] for a given set 
of parameters. In the following, we extend this description 
for the whole range of parameters where such instabilities 
appear. We also discuss in more details than in the 
connections between S and C instabilities, in particular 
through their respective domain of appearance. 

All C behaviors that we observed in the experiments 
have in common to have a large amplitude, and to be 
cyclic, i.e. their trajectory in the phase space follows a 
close cycle. In the time domain, the signal exhibits the 
same pattern, which is repeated indefinitely. However, the 
cadence of the signal is not necessarily periodic, but can 
be erratic. Thus different types of C instabilities may be 
distinguished, and we arbitrarily classified them into three 
groups, that we call Cp, C± and Cs instabilities. 

Among all types of instabilities observed in the MOT, 
Cp instabilities are the most typical deterministic behav- 
ior (Fig. |2| . Indeed, they are characterized by periodic 
oscillations, with a frequency of the order of 1 Hz and a 
large motion amplitude of the order of 100 fim to 1 mm, 
while the population variations are typically 10%. The 
main feature of the cycle is its asymmetry, which can be 
described by the succession of two stages with different 
durations. During the long stage, r and n behave in the 
same way, increasing slowly on a significant amplitude, 
which represents about 30% of the full r amplitude, and 
100% of the full n amplitude. During the short stage, r 
and n change rapidly: n decreases to come back to the 
initial value of the long stage, while r makes a fast oscilla- 
tion, with an amplitude much larger than during the long 
stage. This means that the two stages are not only dif- 
ferent by their duration, but also by the dynamical time 
scales, much faster during the short stage. In fact, the 
characteristic time of the dynamics during the short stage 
is more than one order of magnitude smaller than that in 
the long stage. 

Ci instabilities, illustrated in Fig. 03 corresponds to a 
motion of the cloud very similar to that of Cp instabili- 
ties. Indeed, r exhibits the same behavior along the same 
type of cycle, covered with two different time scales sep- 
arated by one order of magnitude. The main difference 
comes from the erraticity of the motion, which is no more 
periodic: indeed, although the basic pattern of the mo- 
tion remains the same cycle, the duration of each cycle 
fluctuates. From the dynamical point of view, it is con- 
venient use the return time between two cycles, which, in 
chaotic dynamics, is known to be a relevant variable of 
the system^U]. The r return time is constant in periodic 
motions (Cp instabilities), while it is fluctuating in Cj 
instabilities. The n return time is also varying in Ci in- 
stabilities, but other differences as compared to Cp regime 
appear. In particular, the ratio between the fast and slow 
stages of the dynamics fluctuates, so that for some peri- 
ods, the two stages occur with a comparable time scale. 
In short, Cx instabilities appears as Cp instabilities, in 
which noisy fluctuations appear, essentially on the return 
time. 
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Fig. 2. Experimental record of the time evolution of the cloud 
when Cp instabilities occur. In (a), CM motion r; in (b) pop- 
ulation n. Ao = —1.40, h = 10 and I r ep = 1.2 mW/cm 2 . 
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Fig. 3. Experimental record of the time evolution of the cloud 
when Ci instabilities occur. In (a), CM motion r; in (b) pop- 
ulation n. A) = —0.85, h = 10 and I r ep = 1.2 mW/cm 2 . 
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With Cs instabilities (Fig.HJ, any periodicity has dis- 
appeared from the behavior of both r and n. Not only the 
return time is fluctuating, but also the amplitude of the 
cycles changes with time. In fact, speaking of cyclic behav- 
ior in the present case is excessive, except that the basic 
pattern of the motion keeps similarities with the cycle of 
Cp instabilities, in particular the two stages with different 
time scales. However, even the motion cycle is irregular, 
with secondary fast oscillations appearing during the slow 
stage, while the amplitude of the n cycle can fluctuate in 
a ratio of 1:5. 




(a) 



(b) 



Fig. 4. Experimental record of the time evolution of the cloud 
when C's instabilities occur. In (a), CM motion r; in (b) pop- 
ulation n. Aq = —0.75, Ji = 10 and I r ep — 1.2 mW/cm 2 . 



The differences between the three C behaviors are par- 
ticularly well illustrated by the power spectrum of r, as 
shown in Fig. [SJ The spectrum of Cp instabilities exhibits 
a first large and narrow peak, at about 1 Hz, correspond- 
ing to the main period of the signal, followed by a series 
of harmonics (Fig. These harmonics are a signature 
of the second time scale, much faster than the basic pe- 
riod, which appears in the fast oscillation of r. For Cj 
instabilities (Fig. I5t>) , the first peak remains, demonstrat- 
ing that the signal remains essentially periodic. However, 
the regularly spaced harmonics have disappeared, but high 
frequency components remain. They are distributed errat- 
ically, but have globally a larger weight than in Cp insta- 
bilities. Finally, for Cs instabilities, the main frequency 
component has decreased drastically, and the spectrum 
may rather be considered as a wide spectrum, as those 
observed in chaotic or stochastic signals. Unfortunately, it 
is impossible to distinguish between these two possibilities 



through the spectrum analysis of the behavior. To do so, 
it is necessary to turn to more powerful techniques, such 
as the reconstruction of the attractor of the dynamics. 



5-i 



2- 



(a) 




Fig. 5. Spectra of the CM location time evolution in case of 
C instabilities. Figures (a), (b) and (c) correspond respectively 
to Fig. EK, Fig. Hi and Fig. H*. 



Indeed, the MOT is a dissipative system, and any de- 
terministic behavior lies in the phase space on an attrac- 
tor. In the case of deterministic chaos, this attractor is 
complex, but has usually a structured shape, easily rec- 
ognizable. On the contrary, if the behavior is dominated 
by noise, there is no attractor, and the trajectories fill the 
whole phase space. The reconstruction of the attractor of 
a system from experimental time series is a well mastered 
operation. It can be performed following several methods 
(delays, derivatives), and needs usually additional steps, 
as the plot of the Poincare section. In the present case, 
the use of return times is particularly well adapted, as 
it appears as one of the main properties of the behav- 
ior. The return time diagram, which is equivalent to a 
Poincare section, consists in plotting the return time be- 
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tween cycles n and n + 1, as a function of the return time 
between the cycles n — 1 and n^J. The Poincare section 
is a cross-section of the attractor, and has a dimension 
decreased of one unit as compared to the attractor. Thus, 
first return time diagrams of Cp instabilities give just a 
point, as expected from a cyclic behavior. Fig.Elshows the 
first return time diagram in the case of Ci instabilities. It 
is a good illustration of all the return times diagram we 
have obtained for Cj and Cs instabilities: points appear to 
be distributed randomly in the phase space, without any 
structure. Thus we can conclude that the behavior is ei- 
ther high-dimensional deterministic chaos, or a stochastic 
dynamics. From the point of view of the model discussed 
below, these two hypotheses are equivalent, as the model 
includes only three degrees of freedom. 
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tal unstable interval 813 — A3 — A% remains more or less 
constant. When C instabilities merge for 1\ = 41 5, they 
appear on a narrow interval 8 23 > (Fig.©. This inter- 
val increases rapidly until I-y = 7.5/5 an d ^23 = 0.8. For 
I\ > 7.5/5, £23 increases more slowly, to reach the value 
of #23 = 1 in h = 20/s. The value of ^23 depends also 
on the other parameters of the system. Fig.|Hh illustrates, 
as an example, how it depends on the repumper intensity 
Irepi at given I\: 823 varies from for I rep ~ 0.4 mW/cm 2 
to 1 for I rep > 1 mW/cm 2 . 



25 

o 

c 10 



Fig. 6. Typical first return time diagram, obtained from a 
time series of 200 s with Ao = — 1, I\ — 16 and I rep = 7.5 
mW/cm 2 . 



The links between the different C regimes appear clearly 
when one looks at the dependence of the behavior ver- 
sus the different control parameters, and in particular I\ 
and Aq. It was already shown in |S] that for small trap 
beam intensity I\ (typically 1\ < 3), the atomic cloud 
exhibits only S instabilities. When 1\ is increased, S in- 
stabilities still exist, but they are progressively superseded 
by Cinstabilities. The appearance of C instabilities occurs 
progressively, at the cost of S instabilities. For intermedi- 
ate values of 1%, both types of instabilities exist. Their 
typical distribution versus Aq is illustrated in Fig. far 
from resonance, the cloud is stable; as the resonance is ap- 
proached, S instabilities appear for a detuning Aq = A\. 
Then Cinstabilities appear in A2 > A\. If the detuning 
is still increased, C instabilities disappear in A3 at the 
benefit of a stable behavior. Finally, the cloud vanishes in 
A4. As /1 is increased, the interval #23 — A 3 — A 2 where 
C instabilities occur, increases at the cost of the width 
8 \2 = A2 — A\ where S instabilities occur, while the to- 
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Fig. 7. This figure illustrates the evolution of the behavior as a 
function of the detuning for I\ — 6.8 and I rep = 1.5 mW/cm 2 ; 
The full line reports the population, while the dashed lines 
separate the domain of different behaviors: st. stands for stable, 
S (C) for S (C) instabilities. 
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Fig. 8. Width ^23 of the C instabilities area as a function of 
(a) the MOT intensity Ji and (b) the repumper intensity 7 rep . 
In (a), Irep = 1.5 mW/cm 2 , and in (b), I\ — 18 
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Between A2 and A3, the different types of C instabil- 
ities appear, following a constant scenario. In A%, when 
they merge abruply from a stable or S behavior, they are 
Cp instabilities, and as Aq is increased, they transform 
successively in Cj , then Cs instabilities in A 3 . The evo- 
lution is continuous, without abrupt changes, and the Ci 
and Cs instabilities appear rather as two different levels of 
deterioration of Cp instabilities by noise. From this point 
of view, Ci instabilities appear as an intermediate stage 
between Cp and Cs instabilities where noise destroys only 
the periodicity, without affecting the cycle itself. 

The amplitude of the oscillations, of the order of 100 
/im in A2, is much larger than that of the S instabilities 
they merge from, which is typically 30 pm. This is an 
interesting result, because in the most usual bifurcations 
between stable and cyclic behaviors, as the Hopf bifurca- 
tion, the cycle merges progressively from a zero amplitude. 
Such an atypical behavior can be considered as a signa- 
ture of the present system, and must be retrieve in the 
behavior predicted by the model. 

When Ao is increased from A2 , this amplitude still in- 
creases, such that the oscillations reach a 100% contrast in 
z and an amplitude of 40% in n (Fig.0J). Simultaneously, 
the instabilities frequency remains almost constant, as il- 
lustrated on Fig. for different values of the intensity I\. 
The frequency reported here is the main frequency com- 
ponent for Cp and Cx instabilities, and the inverse of a 
mean return time for Cs instabilities. It appears clearly 
that at given I\, the main frequency does not change be- 
tween Ai and A3. This is another interesting result, also 
very untypical in dynamical systems, where the nonlin- 
ear resonance frequencies usually depend stricly on the 
parameters. 

To conclude this section, let us summarize the main 
characteristics of the C instabilities. Three types of be- 
haviors have been identified, depending on their degrees 
of stochasticity: the Cp instabilities are strictly periodic, 
the Ci instabilities remain cyclic, but are no more pe- 
riodic, while finally, the Cs are neither cyclic, nor peri- 
odic. However, this last regime differs drastically from S 
instabilities described in UJ, by their amplitude, and by 
a residue of the cycles they merge from. C behaviors ap- 
pear abruptly, with a non-zero amplitude, and their main 
frequency appears to be independent from the detuning. 



4 Model 

To determine the origin and the exact nature of the in- 
stabilities observed in the experiments, we need to build 
a model able to reproduce also the complex stochastic dy- 
namics observed in [§] . Thus it is logical to use the model 
introduced in [8] and described in details in It is a ID 
model based on the shadow effect induced by the inten- 
sity gradients produced by the absorption of the trapping 
laser beams in the cloud |12U13j . The aim of this model 
is not to reproduce as finely as possible the experimental 
system, but on the contrary to be as simple as possible, 
enlighting the fundamental mechanisms leading to the in- 
stabilities. In its final form, the model reduces to a set of 
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Fig. 9. Experimental evolution of the angular frequency of the 
instabilities as a function of the detuning for different values 
of the beam internsity: in (a), I\ — 4.1; in (b), I\ — 6.8; in (c), 
h = 10; in (d), h = 13.6; in (e), h = 18. I r ep = 1.5 mW/cm 2 
in all cases. 



three autonomous equations, i.e. three equations not de- 
pending explicitely on time. This is important, as it is the 
minimum number of degrees of freedom necessary to gen- 
erate complex dynamics, in particular deterministic chaos. 
The model writes: 



V 



1 



dZ 
~dt 

dV_ 

dt ~ Mv 
dN 
~dt 



■Fi 



= B (1 - Z 2 - N) 



(la) 
(lb) 
(1c) 



where Z = z/zq, V = v/v r and N = n/riQ are the reduced 
variables of the MOT. z and v is the location and the ve- 
locity of the center of mass of the cloud along the unique 
axis z of the system, while n is the number of atoms in- 
side the cloud, zq is a phenomenological size introduced to 
take into account the transverse distribution of the trap 
laser beams, v r is the recoil velocity (y r — hk/m), and rig 
is the equilibrium population of atoms in the cloud. The 
origin of z coincides with the "trap center" , that is, the 
zero of the magnetic field. B is the population relaxation 
rate, M the mass of the cloud and Ft the global force ex- 
erted on the atoms by the two counterpropagating beams. 
To evaluate Fx, we assume a multiple scattering regime, 
i.e. a constant atomic density p in the cloud. Then Ft is 
deduced from the equations of propagation of the beams 
inside the cloud 9 . 

Most of the theoretical parameters are the exact coun- 
terpart of the experimental parameters, as e.g. the mag- 
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netic field gradient or the beam intensities. In this case, we 
used in the model the same values as those of TableQ] It is 
not the case for all parameters, either because of the sim- 
plicity of the model or because they cannot be measured 
easily in the experiment. In particular, no and p cannot 
be accurately evaluated in the experiments. Thus in the 
simulations, they are fixed at experimental averaged val- 
ues, and they have been varied on a wide range to check 
their value is not critical. Finally, to perform the com- 
parison between the experiments and the present model, 
we sometimes need to study the behavior of the system 
when noise is added. This has been done in the same way 
as in [{J], by adding gaussian white noise on I\. Tab. 
summarizes the parameters used in the following. 



Table 2. Parameters used in the numerical simulations. The 
range corresponds to the interval explored numerically, while 
the other sets refer to most of the results presented in this 
paper. 





range 


set #1 


set #2 


G (Gem" 1 ) 


14 


14 


14 


B (s- 1 ) 


3 < B < 30 


3 


3 


h 


2 < Ii < 30 


25 


30 


p (cm -3 ) 


10 10 < p < 3 x 10 10 


2 x 10 10 


2 x 10 10 


S (m 2 ) 


10 -6 < S < 3 x 10~ 6 


10" 6 


10~ 6 


2„ (m) 


10~ 2 < z < 3 x 10 -1 


3 x 10~ 2 


3 x 10~ 2 


no 


10 7 < n < 10 9 


10 8 


6 x 10 s 


A, 


5 < A < -0.5 


-1.5 


-1.5 



5 Stationary solutions 

The model obtained above is described by a set of three 
autonomous equations, and thus could exhibit complex 
behaviors, including periodic and chaotic oscillations, able 
to explain the dynamics observed experimentally. To know 
if such a complex dynamics occurs effectively with our 
parameters, we need first to evaluate the stability of the 
stationary solutions, and thus to calculate the stationary 
solutions themselves. This work has already been partially 
presented in for stable stationary solutions, while we 
are interested here by unstable stationary solutions. How- 
ever, for sake of clarity, we recall here some of the general 
results given in 0, before to start the analysis of the un- 
stable stationary solutions. 

The three stationary solutions Z s , V s and N s are given 
by Eq^ when the left sides are put to zero. As discussed 
m El j V s and N s can be deduced easily from Z s , and thus 
the discussion is reduced to that of Z s , the equation of 
which can be resolved numerically. The global shape of Z s 
is illustrated in Fig. where it is plotted as a function 
of Aq and no- The basic characteristic of this diagram is 
the fold in the stationary solutions, due to several abrupt 
slope changes. The shape of the fold depends on the pa- 
rameters, in particular on no. Fig.llllshows four examples 
corresponding to a situation leading to basically different 



atomic dynamics. For no = 0.5 x 10 8 (fig. Illb). Z s in- 
creases smoothly with Aq (i.e. N s decreases slowly). For 
Ao ~ 0.1 (and thus outside the graph), the cloud vanishes 
through a narrow bistable cycle, where Z s jumps abruptly 
from a value of the order of 0.2 to a value close to 1. As no 
increases, this bistable cycle appears for smaller Z s (and 
thus larger N s ), and becomes physically significant. Fig. 
Illb shows Z s for no = 2.5 x 10 8 and a bistable cycle for 
—0.3 < Ao < —0.25. If no is further increased, the bistable 
cycle disappears, but it remains a fold corresponding to 
two abrupt slope changes of Z s versus Ao (fig. Illb . no 
= 3.4 x 10 8 ). If no is still increased, the fold remains, but 
it becomes smoother (Fig. II Hi for no ~ 4 x 10 s ). 
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Fig. 10. Stationary solutions of equations versus no and 
Ao. The figure represents Z s . Other parameters correspond 
to the set #1 given in table [5] N, F, SN and SF zones (each 
corresponding to different level of greys) describes the nature of 
the fixed point associated with the stationary solution: stable 
Node, stable Focus, Saddle Node and Saddle Focus. 

The results of the linear stability analysis have been 
detailed in [§] . It was shown that the stability and nature 
of the stationary solutions evolve along the fold. In par- 
ticular, the solutions are unstable not only on the central 
branch of the bistable cycle, as it is usual, but also on the 
upper branch of the bistable cycle, and even when there 
is no bistability. This is illustrated in Fig. ^2 where the 
unstable solutions are plotted in a dashed line. As we deal 
here with deterministic instabilities, the interesting situa- 
tion corresponds to Fig. Illb . where the unique stationary 
solution is unstable on the fold. Note the difference with 
the cases studied in where the stationary solution is 
also unique, but stable. In the present case, as no stable 
solution exists, the system exhibits necessarily determin- 
istic instabilities. 

Fig. ^] details the changes in the eigenvalues in this 
situation. Outside the fold (i.e. Ao < A\ or Ao > A4), 
Zs is stable, with one real eigenvalue A r and two complex 
conjugate eigenvalues A ± ito: the fixed point associated 
to the stationary solution in the phase space is a stable 
focus (F zone in Fig. 1101) . The real numbers — A r and —A 
are the damping rates of the stationary solution, and to its 
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Fig. 11. Evolution as a function of the detuning of the sta- 
tionary solution Z a of equations The full (resp. dashed) 
line corresponds to a stable (resp. unstable) solution. In (a), 
n = 0.5 x 10 8 ; in (b) n = 2.5 x 10 8 ; in (c), n = 3.4 x 10 8 ; in 
(d), no — 4 x 10 s Other parameters correspond to the set #1 
of Table H 



eigenfrequency. The transition from stable focus solution 
to unstable saddle focus solution occurs in A\ and A±, 
through a Hopf bifurcation, where A = and w ^ 0. As it 
is usual for such a bifurcation, we expect to observe, on the 
unstable side, a cloud moving on a limit cycle with a fre- 
quency u>, i.e. in the present case of the order of 200 s _1 (30 
Hz) for both Hopf bifurcations. Another transition, from 
a saddle focus solution to a saddle node solution, occurs 
in A2 and A^, when to vanishes. The behavior here cannot 
be deduced from the stationary solutions, and numerical 
simulations are necessary. Deterministic instabilities are 
expected to occur between Ai and A4. In all other areas, 
the stationary solutions are stable, and therefore, deter- 
ministic instabilities cannot occur|5]. In the next section, 
we detail the results of numerical simulations performed 
in the unstable area. 



However, before to discuss in detail the dynamical be- 
haviors predicted by the model in the different situations, 
let us look at the influence of the other parameters on the 
stationary solutions. As discussed in the previous section, 
we must distinguish between theoretical parameters with 
an exact experimental counterpart, as I±, for which the 
comparison with experiments is direct, from those without 
an exact experimental counterpart, for which the analysis 
is more delicate. It is in particular the case for no and 
p, which are both linked to I rep and the vapour pressure. 
Finally, the influence of B and z$ should be checked, as 
their experimental determination is unprecise. 

Fig. El illustrates the role of the I\ value on the be- 
havior of the cloud: it represents the stationary solution 
Z$ versus the detuning, for different values of the inten- 
sity. The figure shows that I± acts as no on Zs'. an in- 
crease of Ii makes the fold steeper, and eventually leads 
to bistability. However, some more subtle changes occur, 
as illustrated by Fig. O where Zs have been plotted as 
a function of Aq for different values of no, as in Fig. 1111 
but for a smaller intensity. In these new conditions, the 
intermediate area between the stable fold and bistability, 
where the stationary solution is unique and unstable, has 
almost disappear. This result may be generalized: in the 
simulations, we observed that the area corresponding to 
an unstable unique solution disappears for small intensi- 
ties (typically I\ < 10). This result is in agreement with 
the experimental observation that C instabilities exist only 
for large intensities. 
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Fig. 13. Evolution as a function of the detuning of the sta- 
tionary solution Z a of equations for different values of the 
intensity I\. The full (resp. dashed) line corresponds to a sta- 
ble (resp. unstable) solution. In (a), / = 10; in (b) / = 15; 
in (c), I = 20; in (d), I = 25; in (e), / = 30; in (d), / = 35. 
Other parameters correspond to the set #1 of Table |21 with 
n = 2 x 10 s . 



Fig. 12. Evolution as a function of the detuning of the sta- 
tionary solution Zs and its eigenvalues, for the parameters of 
Fig. lilt . The stationary solution is given through the full (sta- 
ble) and dashed (unstable) bold lines. The dashed line noted 
uj represents the imaginary part of the complex eigenvalues, 
while the full lines correspond to the real part of the three 
eigenvalues. The full line remaining negative corresponds to 
the eigenvalue which is real everywhere. 



The atomic density value used in the present simula- 
tion is p — 2 x 10 10 cm~ 3 , which represents an average of 
the density we measured experimentally. To evaluate the 
influence of this value on the predicted behavior, we have 
plotted in Fig. ^| the evolution of Zs for different values 
of p. As for no and I± , the Zg curve evolves from an almost 
flat dependence for large p, towards bistability for small p, 
with an intermediate SF zone. This can appear as surpris- 
ing, because it seems to mean that nonlinear behaviors, 
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Fig. 14. Evolution as a function of the detuning of the sta- 
tionary solution Z a of equations Q for different values of no 
and a smaller intensity than in Fig. 1111 The full (resp. dashed) 
line corresponds to a stable (resp. unstable) solution. In (a), 
no = 0.25 x 10 8 ; in (b) n = 0.5 x 10 s ; in (c), n = 0.7 x 10 s ; 
in (d), no = 1 x 10 8 ; in (e), no = 2 x 10 s . Other parameters 
correspond to the set #1 of Table |5] with Ji = 15. 



corresponding to the bistable cycle, need a small atomic 
density. In fact, this reasoning is false, because it does not 
take into account the role of the other parameters, in par- 
ticular uq, which is able to compensate for the variation 
of p. For example, fig. H6l shows the (Aq, no) diagram for a 
smaller p value than in Fig. 1111 it has the same properties 
as that in fig. 1111 except that the population are much 
larger. However, one remarks slight differences, in partic- 
ular a small decreasing of the Z values, and first of all 
a small decreasing of the SF zone width. This is another 
general result: in the simulations, when the atomic den- 
sity is decreased, the Z s curves globally flatten, so that 
the fold becomes less steep. Thus all instabilities disap- 
pear for very small densities. This is in agreement with 
the experimental results illustrated in Fig. |SJd, where the 
unstable interval width is reported as a function of I rep : 
the instabilities disappear for small I rep , i.e. when the ef- 
ficiency of the repumper - and thus the atomic density - 
decreases. 

As detailed in [2], the value of zq = 3 cm used in the 
experiments has been evaluated from the trap beam waist, 
taking into account the beam intensity as compared to 
the saturation intensity. As for p, we want to evaluate 
how critical is this choice by plotting Zs versus Aq for 
different values of z (Fig. El- It appears clearly that a 
decreasing of zo corresponds to a shift of the fold and 
the bistable cycle towards resonance. Thus, for zo values 
smaller than 3 cm, the discrepancy between simulations 
and experiments increases quantitatively, as instabilities 
will appear at smaller detuning. For values smaller than 1 
cm, the change is drastic, as the SF and bistable zones, and 
thus instabilities, disappear. On the contrary, for larger zq, 
the bistable cycle widens and shifts off resonance. 

The last parameter to be considered is B. This param- 
eter appears to be not critical at all, and in particular, a 
change from e.g. B = 5 s" 1 to B — 1 s" 1 does not change 
in a noticeable way the values of Zs- The main change 
concerns the smaller real eigenvalue, which takes typically 
the value of —B. However, this eigenvalue plays a minor 
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Fig. 15. Evolution as a function of the detuning of the sta- 
tionary solution Z a of equations Q for different values of the 
atomic density. The full (resp. dashed) line corresponds to a 
stable (resp. unstable) solution. In (a), p = 1 X 10 10 cm" 3 ; in 
(b), p = 1.5 X 10 10 cm" 3 ; in (c), p = 2 x 10 10 cm" 3 ; in (d), 
p = 2.5 x 10 10 cm" 3 ; in (e), p = 3 x 10 10 cm" 3 . Other param- 
eters correspond to the set #1 of Table with no = 2 x 10 8 
and I\ — 25. 
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Fig. 16. Evolution as a function of the detuning of the station- 
ary solution Z 3 of equations for different values of no and a 
smaller atomic density than in Fie. till The full (resp. dashed) 
line corresponds to a stable (resp. unstable) solution. In (a), 
n = 2 x 10 s ; in (b) n = 4 x 10 s ; in (c), n = 6 x 10 8 ; in (d), 
n = 8x 10 8 ; in (e), n = 1 x 10 9 ; in (f), n = 1.2 x 10 9 ; in (g), 
no = 1.4 x 10 9 . Other parameters correspond to the set #1 of 
Table with p = 1 x 10 10 cm" 3 . 

role in the dynamics, as it remains always real negative, 
and thus this change has a negligible effect on the dynam- 
ics. 

It appears from the above analysis that the existence 
of the unstable area does not depend critically on the val- 
ues of no, 1%, p, zq and B. In particular, the relative poor 
accuracy in the knowledge of the experimental values of 
some parameters, as no, p or Zq, does not appear as a limi- 
tation in the above study, because a change of some tens of 
percents around the default values used in the simulations 
does not alter the results. The numerous approximations 
at the origin of the model lead probably to larger errors. 

6 The unstable fold: deterministic instabilities 

As shown in the previous section, it exists a range of pa- 
rameters where the stationary solution is unique and un- 
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Fig. 17. Evolution as a function of the detuning of the station- 
ary solution Z 3 of equations for different values of zo. The 
full (resp. dashed) line corresponds to a stable (resp. unstable) 
solution. In (a), Zo = 1 cm; in (b), zo = 1.1 cm; in (c), zo = 1.5 
cm; in (d), zo = 2 cm; in (e), zq — 3 cm. Other parameters 
correspond to the set #1 of Table [5] with no = 2 x 10 8 and 
h = 25. 



stable. Such a situation leads inevitably to determinis- 
tic instabilities, with shape and characteristics obtained 
through numerical simulations of the model for the cor- 
responding sets of parameters. In the present section, we 
discuss the behavior obtained by such simulations, in a 
situation similar to Figs II lb and 1 121 but for a slightly dif- 
ferent set of parameters (set #2 of Tab.[2J. Fig. 1181 shows 
for these conditions the evolution of the eigenvalues as 
a function of the detuning. The sequence is identical to 
that followed in Fig. but the unstable zone is wider. 
Starting off resonance, a Hopf bifurcation occurs in A\. 
In Z\ 2 , the eigenvalues become real, and thus the eigen- 
frcquency disappears, until ^3, where two eigenvalues are 
again complex. Finally, a second Hopf bifurcation occurs 
in A4, and the stationary solutions become stable again. 
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Fig. 18. Evolution as a function of the detuning of the eigen- 
values of the stationary solutions. The full line, which is a 
plot of the real part A of the eigenvalues when they are pos- 
itive, put in evidence two bifurcations, in Ai = —0.4024 and 
A4 = —0.3566. The dashed line represents the corresponding 
imaginary part u (i.e. the eigenfrequencies). The eigenvalues 
are real between A2 = —0.3981 and A3 = —0.3760. Parameters 
corresponds to the set #2 of Tab. [2] 



6.1 The behavior in the vicinity of the Hopf bifurcation 

To understand the origin of the C instabilities, we have 
studied in detail the behavior of the model in the close 
vicinity of the Hi Hopf bifurcation, by varying Aq slightly 
above A\. The expected scenario has been extensively de- 
scribed in the litterature^UI- O n the left of Hi, the sta- 
tionary solution is stable, and thus the behavior is sta- 
tionary. In Hi, the solution becomes unstable, but a limit 
cycle merges from the unstable fixed point: the behav- 
ior becomes periodic, with a zero amplitude in Hi, and 
a frequency corresponding to the relaxation frequency of 
the unstable solution. On the right of the bifurcation, the 
amplitude of the oscillations grows, while the frequency 
remains the same as the eigenfrequency in the vicinity 
of Hi. Usually, when the control parameter is increased 
from Hi , the cycle shape changes progressively, while the 
interval between the oscillation and relaxation frequencies 
becomes larger. This standard scenario is absolutely not 
followed by the present model. On the contrary, several 
abrupt changes in the behavior occur in a very narrow in- 
terval of Aq, leading from a classic regular limit cycle to 
the Cp instabilities. 

In Hi appears, as expected, a limit cycle. Figs. HHl and 
1201 show the evolution of this limit cycle on the unstable 
side of the bifurcation between A\ and A\ — —0.4010. 
Note that the explored interval is so narrow (1.3 x 10~ 3 ) 
that an experimental observation of the described phe- 
nomena cannot be considered. Following the Z coordinate, 
the cycle amplitude grows rapidly to reach a value of typ- 
ically 10% of Zs, while the amplitude on N remains very 
small (0.1% of Ns). The cycle remains relatively well cen- 
tered on Zs, but is shifted compared to Ns, such that Ns 
is well outside the cycle. This is not an exceptional situa- 
tion, as it simply means that the basin of attraction of the 
cycle is curved in the vicinity of the unstable fixed point. 
Another characteristics of the limit cycle is its frequency, 
which remains of the order of magnitude of the eigenfre- 
quency. For example, for Aq — —0.4015 (Fig. 119b). the 
behavior frequency is 30 Hz, for an eigenvalue of 27 Hz. 
Thus the global behavior in the interval (Ai, A^) appears 
to be the usual one in the vicinity of a Hopf bifurcation. 

However, for A > A[, the limit cycle becomes unsta- 
ble and is replaced by another periodic orbit, with a much 
more complex shape (Fig. 12 II) and a much longer period. 
The amplitude is almost 5 times larger for Z and more 
than 10 times for N. The trajectory consists in several dif- 
ferent stages: a diverging spiral off the fixed point, followed 
by a large loop and a convergent spiral until the fixed 
point. The frequencies of the two oscillating stages are dif- 
ferent: this is not surprising, as the diverging one is clearly 
linked to the fixed point, and thus to its eigenfrequency, 
contrary to the convergent one. One finds effectively a fre- 
quency of 24.4 Hz for the diverging spiral, corresponding 
exactly to the eigenfrequency (oj = 2ir x 23.5 Hz), while the 
frequency of the converging spiral is 77 Hz. However, the 
main frequency of the behavior is 2.6 Hz, i.e. one order of 
magnitude slower than that of the Hopf cycle. The prop- 
erties of the trajectories, in particular the tangency to the 
unstable point and the large loop in the phase space, are 
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Fig. 19. Time evolution of the atomic cloud location Z in 
the close vicinity of the A\ bifurcation, on the unstable side. 
The dashed line represents the value of the unstable stationary 
solution Z s . In (a), A = -0.4023; in (b), A = -0.4020; in 
(c), A = -0.4015; in (d), A = -0.4010. Other parameters 
are those of set #2 of Tab.H 
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Fig. 20. Time evolution of the atomic cloud population N in 
the close vicinity of the A\ bifurcation, on the unstable side. 
The dashed line represents the value of the unstable stationary 
solution Ns- Parameters are the same as in Fig 1191 



characteristic from a homoclinic behavior, when the sta- 
ble and unstable manifolds of the fixed point are almost 
connected. A accurate analysis of these manifolds would 
be necessary to conclude about this point. Note that in 
a very narrow interval around A^, generalized bistability 
occurs between the Hopf cycle and the homoclinic one. 

When Aq is increased from A[, the cloud exhibits pe- 
riod doubling and chaos (Fig. 122(1 . The trajectories keep 
the same shape, in particular with the two spiraling episodes 
and the large loop, but the periodicity is modified or is 
lost. For example, when the period is double, variations 
appear essentially on the amplitude of the loop together 
with that of the diverging oscillations (Fig. 122b). In the 
chaotic zone, the irregularities appear also on these am- 
plitudes, but bursting events appear sometimes between 
these two stages (Fig. We did not perform a precise 
analysis of these behaviors, mainly because they appear 
on a so narrow interval that there is no chance to ob- 
serve them experimentally. Indeed, chaos disappears for 
A > A'l, with A'l = -0.4005, and thus the homoclinic 
behavior appears on an interval of 5 x 10~ 4 . However, a 
simple test can be done by reconstructing the attractor 
of the dynamics fFig. I23fl . A glance at the result shows a 
definite structure, and not random distributed points, and 




Fig. 21. Time evolution of the atomic cloud location Z in (a), 
and population N in (b), for Ao = —0.4010. Other parameters 
are the same as set #2 in Tab. [5] 
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thus confirms that this behavior presents all the charac- 
teristics of deterministic chaos. 
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Fig. 22. Time evolution of the atomic cloud location Z in 
the vicinity of the chaotic zone. In (a), for Ao — —0.4010, the 
motion is periodic; in (b), for Ao = —0.4006, period doubling 
appears; in (c), for Ao = —0.4005, the motion is chaotic. Other 
parameters are the same as set #2 in Tab.|5] 



6.2 Cp instabilities 

For Ao > A'{, the homoclinic behavior disappears, and a 
new type of periodic instabilities appear (Fig. I24fl . There 
is no fundamental difference between the homoclinic in- 
stabilities and the present behavior, except that the latter 
has a physical meaning, as it appears in the simulations on 
a significant Z\ interval, about 0.04.T for the present pa- 
rameters. One observes in Fig. 124b , the same three stages 
as in Fig. 1221 which means that the origin of this behav- 
ior is the same as for the homoclinic instabilities. How- 
ever, these three stages exist only in the close vicinity of 
A'{: when Aq is increased, the diverging helix around the 
fixed point disappears, and only the two stages indepen- 
dent from the fixed point remain (Fig. 124b and c). This 
means that in this new behavior, the trajectories never 
approach the fixed point, and thus the dynamics does not 
depend on the local properties of the fixed point. 

Each period of the new behavior may be divided in two 
stages with different durations. During the fast stage, Z 
makes an oscillation, first growing then decreasing, while 
N decreases; during the slow stage, Z and N grow. This 
is similar to that observed in experiments with Cp in- 
stabilities, and thus it is interesting to check if the other 
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Fig. 23. Chaotic trajectories of the cloud motion in the (Z, N) 
phase space. Ao = —0.4005. Other parameters are the same as 
in Tab.0 
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Fig. 24. Time evolution of the atomic cloud location Z 
(left) and population N (right) in the unstable zone. In (a), 
Ao = -0.4004; in (b), Ao = -0.39; in (c), A Q = -0.37. Other 
parameters are the same as in Tab. [5] 



properties of Cp instabilities can be retrieve in the present 
dynamics. 

Two untypical properties were noticed in the Cp in- 
stabilities, concerning the non zero amplitude of the os- 
cillations when they appear, and their almost constant 
frequency along their interval of existence. The amplitude 
of the oscillations as a function of Aq is plotted in fig. |2*51 
In Aq = —0.4, when the instabilities appear, their ampli- 
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tude is already almost 0.3. In fact, because the Hopf limit 
cycle exists on a very narrow interval, this is not strictly 
true. But from a physical point of view, it is clear that 
instabilities appear with a non zero amplitude, as in the 
experiments. Concerning the frequency, its value tends to 
zero in the vicinity of A", but it increases rapidly to reach 
a value of the order of 10 Hz, and remains between 10 Hz 
and 13 Hz on most of the unstable interval, as shown in 
Fig. 1261 fthe behavior for Aq > —0.36 is discussed below). 
Finally, to complete the comparison with the experiments, 
Fig. 1271 shows the power spectra of Z. In the vicinity of 
Hi, it is very characteristic, with a main frequency and 
large amplitude harmonics decreasing progressively (Fig. 
127b ) , as in the experiments (Fig. [SJ) . As the detuning is 
increased, the amplitude of the harmonics decreases, and 
several new frequencies appear in the spectrum, but each 
component has individually an amplitude negligible com- 
pared to the main frequency (Fig. EZt> and c). 
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Fig. 25. Amplitude of the dynamics versus detuning. The 
squares, discs and triangles corresponds respectively to the 
oscillation amplitude, the maximum value and the minimum 
value reached by z. Parameters are those of Tab. [5] 
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Fig. 26. Evolution as a function of the detuning of the insta- 
bilities frequency. 



To summarize, we found that the periodic instabili- 
ties in the vicinity of Hi have the same shape, spectrum, 
amplitude evolution and frequency evolution as the Cp 
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Fig. 27. Power spectra of the atomic cloud location Z in the 
unstable zone. Parameters are the same as in Fig l24l 



instabilities observed in the experiments; all these points 
confirm that the model reproduces here the Cp instabili- 
ties (Fig. Hand OH). 

However, several quantitative discrepancies appear be- 
tween the present results and the experimental observa- 
tions, as e.g. for the Aq values or the instabilities fre- 
quency. These differences are relatively small, except for 
the Aq interval where Cp instabilities exist: it is typically 
of the order of 1 in the experiments, while it is smaller than 
0.1 in the simulations. This last value could be increased 
by increasing the value of I\ in the simulations, reducing 
the difference to less than a factor 10. Considering the 
extreme simplicity of our model and its numerous approx- 
imations, it is clear that it is able to reproduce strikingly 
the Cp instabilities, with a surprisingly good agreement. 



6.3 C'i instabilities 

In the experiments, Cp instabilities are replaced, as the 
detuning is changed, by Cj instabilities. The difference be- 
tween the two regimes appears mainly on the n behavior, 
and in particular in the shape of the signal, which loses its 
regularity. The Ci behavior is not observed in the present 
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model. However, as discussed in the experimental section, 
the Cj instabilities do not appear as a new deterministic 
regime, but rather as Cp instabilities slightly altered by 
noise. Thus, to test the ability of the present model to re- 
produce this behavior, it is necessary to add noise in the 
model. The results obtained in this case are discussed in 
the next section. 



6.4 C,s instabilities 

Fig. |23 shows that the amplitude of the oscillations in- 
creases with Ao, so that the maximum value explored by 
z becomes larger and larger as Aq is increased. As a conse- 
quence, the maximum value reached by z also increases, so 
that finally, the most distant atoms from the trap center, 
situated in z + Az/2, where Az is the size of the cloud, 
reach the border of the trap, in Zq. In the model, these 
atoms are considered to be lost, and thus are subtracted 
to the total number of atoms in the trap. Therefore, a 
new process with a zero characteristic time appears in the 
model through this instantaneous decreasing of n. This 
new process leads to an immediate change of the dynam- 
ics frequency. This is illustrated in fig. 1261 where the tran- 
sition occurs in A = —0.361. For these parameters, the 
frequency is divided by more than a factor 3, decreasing 
from 12 Hz to 3.7 Hz. The new dynamics is illustrated in 
Fig- EH Although the global shape seems to be similar to 
the previous Cp instabilities, a drastic difference appears 
on the dynamics of N, in particular concerning its oscil- 
lation amplitude. While the variations of N as a function 
of time were small in the Cp regime, they appear to be 
much larger in the present regime: in fig. 128b . N varies 
on half of the total interval of values that it can take, 
and when Aq is still increased, this ratio reach 80%, with 
an oscillation from to 0.8 (Fig. 128b). In the latter, the 
cloud empties completely every period, and then fill up 
progressively. This explain the large period of the regime, 
due to a longer time necessary to fill the cloud. This ap- 
pears clearly when fig. El and fig. EEI & re compared: the 
increasing of the period does not correspond to a global 
stretching of the dynamics, as the large oscillation of Z re- 
mains on the same time scale, but rather is a consequence 
of the increasing of the interval between two oscillations. 

The spectrum confirms that the main frequency has 
decreased (Fig. l29|) . However, this is coupled with the ap- 
pearance of large amplitude harmonics, decreasing slowly, 
so that components with higher frequency than in the Cp 
regime keep a significant weight. 

This behavior has several common points with the Cs 
instabilities described in the experimental section. In both 
cases, the regime is the continuation of the Cp instabilities 
when the resonance is approached, the amplitude of the z 
oscillation has a 100% contrast, that of n are much larger 
than in Cp regime, the shapes are similar, and new fre- 
quencies of higher value appear. But the regime obtained 
in the simulations remains periodic, contrary to that ob- 
served experimentally. However, as discussed in the exper- 
imental section, the origin of the crraticity observed in the 
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Fig. 28. Time evolution of the atomic cloud location Z (left) 
and population N (right) in the unstable zone. In (a), Ao = 
-0.36; in (b), A = -0.35. 
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Fig. 29. Power spectra of the atomic cloud location Z in the 
unstable zone. Parameters are the same as in Fie l24l 



experiments could be stochastic, rather than determinis- 
tic. Thus we can hope that the addition on noise in the 
model will transform the dynamics to reproduce the ex- 
perimental results. This influence of noise on the dynamics 
is studied in the next section. 
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7 The effect of noise 

Noise is known to be able to alter drastically the deter- 
ministic behavior of the MOT: in Sj, it has been shown 
that a stationary behavior may be transformed, under the 
influence of noise, in a behavior similar to instabilities. 
Thus a complete study of the behaviors predicted by the 
present model must consider the possible alterations in- 
duced by noise on the dynamics. We present successively 
in this section the results obtained on Cp and Cs insta- 
bilities. 

Fig. |201 illustrates the effect of noise on Cp instabili- 
ties, through the example of the regime plotted on Fig.lM: 
with 2% of noise on I\. Although the behavior becomes 
less regular, with fluctuating amplitudes and a ruffling of 
the small secondary oscillations, the global behavior is un- 
changed, with a still rather regular period and the same 
global shape as without noise. Thus Cp instabilities ap- 
pear to be robust against noise. However, if the response 
to noise is studied in more details, a slight increase of the 
noise influence appears when the Cs area is approached. 
This global behavior allows us to interprete both the Cp 
and the Cj instabilities, which appear in fact to be the 
same dynamics, affected differently by noise: between A2 
and A3, i.e. far from the Cs area, the Cp instabilities are 
very robust against noise, and the presence of technical 
noise in the experiment does not alter neither their shape 
nor their periodicity. As the Cs area is approached, the 
sensitivity of the Cp behavior to noise slightly increases, 
and, although the main characteristics of the Cp instabil- 
ities remain unchanged, the shape and periodicity of the 
regime are affected enough to give the feeling of a new 
regime, namely Cj instabilities. Thus Ci instabilities ap- 
pear, as already suspected in the experimental section, as 
a Cp regime perturbed by noise 
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Fig. 30. Time evolution of Z for the same parameters as in 
Fig. 124b . but when noise is added on I\, The noise level is 2%. 



Fig. shows how noise alters the dynamics of Cs in- 
stabilities. Although the amount of noise is the same as 
for Cp instabilities, it is clear that here, the dynamics 
is deeply transformed. Concerning the Z dynamics, the 
shape of the oscillations remains almost unchanged, but 
the periodicity is drastically altered: indeed, the return 
time of the pulses varies randomly on a range larger than 
100%. The explanation of these large fluctuations on the 



return time of the Z pulses comes from the N dynam- 
ics: here, the fluctuations appear on the amplitude of the 
variable. As the return time is connected to the recon- 
struction time of the cloud, it is logical that fluctuations 
in the initial population lead to fluctuations in the return 
time of Z. The strength of the effect is due, as for the 
stochastic dynamics described in [H] , to an effect of ampli- 
fication of noise, but through a different mechanism than 
that described in Indeed, in the vicinity of zq, noise 
modifies the state of the cloud just before the brutal de- 
creasing of the population: it will be able to slow down 
the crossing of zq, or on the contrary to quicken it. The 
consequence on the number of atoms lost in the process is 
immediate, leading to the large fluctuations of N that can 
be seen in Fig. 1311 If the resulting dynamics is compared 
with the experimental one illustrated in Fig.^ the similar- 
ities between both dynamics appear clearly: the common 
points discussed in the previous section remain, and the 
discrepancies disappear. In particular, the fluctuations in 
the return time of the z pulses, together with those in 
the amplitude of n, are now present: it is clear that we 
reproduce here the Cs instabilities. 
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Fig. 31. Time evolution of the (a) atomic cloud location and 
(b) population, for the same parameters as in Fig. 1281 but when 
noise is added on I\. The noise level is 2%. 



This concludes the present study on the deterministic 
instabilities of the MOT, as we have been able to retrieve 
with our model all the behaviors observed experimentally. 
In particular, the dynamics that appeared as erratic is 
shown to be a deterministic periodic behavior perturbed 
by noise. This allows us to do the link with the studies 
reported in [5] , as we show in this section that noise plays 
once again a key role in the dynamics of the MOT, al- 
though here the basis of the dynamics is deterministic. 



16 Andrea di Stefano, Philippe Verkerk and Daniel Hennequin: Deterministic instabilities in the magneto-optical trap 



8 Conclusion 

The behavior of the cloud of cold atoms produced by a 
magneto-optical trap exhibits a rich variety of dynamics, 
that are well reproduce by a simple ID model described by 
a set of three autonomous ordinary differential equations. 
In |S] were described a set of noise- induced instabilities, 
linked to the topology of the stationary solutions. Here we 
show that experimentally, three different regimes of deter- 
ministic instabilities may also be observed, depending on 
the parameters of the MOT. Some of these regimes appear 
to be purely deterministic (Cp instabilities), some appear 
to be a mixture of deterministic instabilities and effects 
of noise (Ci and Cs instabilities). Theoretically, the same 
model as in 9 allows us to show that the existence of 
these deterministic instabilities are a direct consequence 
of the topological properties that induce for other param- 
eters the stochastic instabilities studied in |U] . This model 
shows that C± instabilities are just Cp instabilities per- 
turbed by noise, while Cs instabilities are a nother de- 
terministic regime appearing when the border of the trap 
beams are reached. This last regime is particularly sen- 
sitive to noise, and the resulting behavior appears as a 
deterministic instability highly perturbed by noise. Thus 
the present study confirms that noise plays a crucial role 
in the dynamics of the atomic cloud. 

The present results have been obtained with a very 
simple model. Although the agreement with the experi- 
ments is surprisingly good, it is difficult to make qunati- 
tative comparisons between such a ID model and a 3D 
experiment. Therefore the next theoretical step would be 
to develop a 3D-model, where some of the parameters of 
the present model, as e.g. the cloud volume, become a 
function of the dynamics variables. 

An interesting perspective of the present results is to 
study the possibility to take advantage of the existence 
of deterministic instabilities in the MOT. In particular, 
if a set of parameters could be found to widen enough 
the chaotic area, the techniques of control of chaos could 
be apply to reach various states that are not accessible 
otherwise, as e.g. denser or colder states. But even periodic 
behaviors can give new interesting informations about the 
MOT physics. Indeed, a complex dynamics covers a larger 
part of its phase space, and in return makes possible the 
determination of parameter values masked in stationary 
behaviors. More generally, a complex behavior enables the 
access to more information about its system, and appears 
usually as a good starting point for a better understanding 
of it. 
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